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ABSTRACT 

The clumping of massive star winds is an established paradigm confirmed by multiple lines of 
evidence and supported by stellar wind theory. The purpose of this paper is to bridge the gap 
between detailed models of inhomogeneous stellar winds in single stars and the phenomeno- 
logical description of donor winds in supergiant high-mass X-ray binaries (HMXBs). We use 
results from time-dependent hydrodynamical models of the instability in the line-driven wind 
of a massive supergiant star to derive the time-dependent accretion rate onto a compact object 
in the Bondi-Hoyle-Lyttleton approximation. The strong density and velocity fluctuations in 
the wind result in strong variability of the synthetic X-ray light curves. The model predicts a 
large scale X-ray variability, up to eight orders of magnitude, on relatively short timescales. 
The apparent lack of evidence for such strong variability in the observed HMXBs indicates 
that the details of accretion process act to reduce the variability due to the stellar wind velocity 
and density jumps. 

We study the absorption of X-rays in the clumped stellar wind by means of a 2-D stochas- 
tic wind model. The monochromatic absorption in cool stellar wind in dependence on orbital 
phase is computed for realistic stellar wind opacity. We find that absorption of X-rays changes 
strongly at different orbital phases. The degree of the variability due to the absorption in the 
wind depends on the shape of the wind clumps and is stronger in case of oblate clumps. 

We address the photoionization in the clumped wind, and show that the degree of ioniza- 
tion is affected by the wind clumping. A correction factor for the photoionization parameter 
is derived. It is shown that the photoionization parameter is reduced by a factor X compared 
to the smooth wind models with the same mass-loss rate, where X is the wind inhomogeneity 
parameter We conclude that wind clumping must also be taken into account when comparing 
the observed and model spectra of the photoionized stellar wind. 
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1 INTRODUCTION 

Massive luminous OB-type stars possess strong stellar winds. The 
winds are fast, with typical velocities up to 2500 km s^ ^ , and dense, 
with mass-loss rates M>10~^ Mq yr^^. The driving mechanism 
for the mass-loss from OB st ars has been ident ified with radia- 
tion pressure on spectral lines dCastor et al.lll975l , CAK). Early on 
lLucv&Solomonl l ll970h suggested that the stationary solution for 
a line-driven wind is unstable leading to shocks in stellar winds. 
lLucv&Whit3 ( ll980l) pr oposed that the winds break up into a pop- 
ulation of dense blobs. lOwocki et al. I l ll988l) performed numeri- 
cal simulations of the nonlinear evolution of wind instabilities and 
showed that even small-amplitude perturbations in the radial ve- 
locity near the wind base result in a highly structured wind, with 
relatively slow, dense shells separated from regions of high-speed, 
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rarefied flow by strong shocks. [peldmeier et al.l Il997al lbh used hy- 
drodynamic simulations to model the evolution of wind instabilities 
and predicted the X-ray emission from single stars matching the 
observed X-ray fluxes. These non-stationary hydrodynam ic simu- 
lations are used in this paper iRunacres & Owockil ( |2002|) studied 
the evolution of wind structures and demonstrated that the winds 
can remain inhomogeneo us out to distances of 1000 R*, while 
iDessart & Owockil ( l2005h presented first attempts of 2-D radiative 
hydrodynamic models of stellar winds. 

These theoretical predictions of stellar wind clumping are con- 
fir med by obser vations. From analysis of wind lines in the UV spec- 
tra |Luc3 jl982h postulated the existence of non-monotonic veloc- 
ity fields in the winds. Hillier ( 199ll) induced the existence of wind 
clumping from the shape of electron-scattering wings in emission 
lines of Wolf-Rayet (WR) stars. Stochastic variable structures in 
th e Hell A4686A emiss ion line in an O-supergiant were found 
bv lEversberg et al] ( Il998l) . and explained as excess emission from 
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the wind clumps. iMarkova et alj l l20()^ investigated the line-profile 
variability of Ha for a large sample of 0-type supergiants. They 
concluded that the observed variability can be ex plained by a wind 
model consi sting of coherent or broken shells. iLepine & Moffall 
( I1999II2OO8I) presented d irect spectroscopic evid ence of clumping 
in O and WR star winds. IPrinia & Ma"ss3 bOld) established spec- 
troscopic signatures of the wide-spread existence of wind clumping 
in B supergiants. 

The inhom o geneit y affects stellar wind diagnostics. 
IStewart & FabianI ( Il98lh used Einstein spectra of 0-stars to 
study the transfer of X-rays through a uniform stellar wind. By 
matching the model and the observed X-ray spectrum of the 
O-type supergiant C, Pup they found that the mass-loss rate derived 
from X-ray spectra is lower by factor of a few than the mass-loss 
rates obtained from fitting the Ha line as well as from observations 
at radio and IR wavelengths. As most plausible explanation for 
this discrepa ncy they suggested the negl ect of the clumping in 
O-star winds. IWaldron & Cassinellil ( l200lh proposed that a clumpy 
wind structure may help to explain the observed X-ray emission 
line profiles in the spectra of O-supergiants. A theory of X-ray 
transfer in clumped wins that accounts for clumps of any shapes 
and optical dep t hs wa s first developed by Feldmeier et al. (2003); 
lOskinovaetail ( |2004|) . Indeed, accounting for wind clumping 
allows to c onsistently model the UV/optical and the X-ray spectra 
of O-stars ( IOskinovaetalj2006h . 

Presently, clumping in first approximation ("microclump- 
ing") is taken into account in all up-to-date stellar wind codes. 
This approximation assumes that clumps are optically thin and 
that the interclump medium is void. Microclumping leads to 
smaller empirical mass-loss rates when using diagnostics that 
depend on processes scaling with density squared (p^), such 
as Ha and other emission lines, or radio free-free emis sion 
(lHillieJl99ll;lHamann & Koesterkell9"93) . lBouret etai] ( l2005h and 
iFullerton et all ( l2006h demonstrated that mass-loss rates must be 
reduced by orders of magnitude in order to reproduce the UV- 
resonant lines (which depend on density linearly) simultaneously 
with the based diagnostics^. These drastic reductions of empir- 
ically deduced mass-loss rates would have strong consequences for 

the stellar evolution and stellar feedback models. 

A way to solve this problem was suggested bv lOskinova et al.l 

( I2OO7I) . It was found that waiving the microclumping approxima- 
tion i.e. allowing for clumps of arbitrary optical depth ("macro- 
clumping") in the spectral models eliminates the need to strongly 
reduce stellar mass-loss rates. It was shown that the strength of 
optically thick diagnostic lines, such as the P V resonance line in 
OB star spectra, is reduced by macroclumping effects, while the 
strength of optically thin lines, such as Ha, is not affected. There- 
fore, accounting for macroclumping in analyzing stellar spectra 
yields empirically derived mass-loss rates that are only a factor 2 - 
3 lower than the "classical" CAK mass-lo ss rates. These conclu- 
sions are confirmed by 2D model studies l lSundgvist e7ai]|2010l) 
and full 3D radiative-transfer models (Sudan et al., submitted). 

The CAK mass-loss rates of massive stars are in good agree- 
ment with observed X-ray fluxes from high-mass X-ray binaries 
(HMXBs). HMXBs are the produc ts of binary star evolution (e.g. 
IShklovskvl Il967l : llben et al.l [T99I ). These systems consist of an 
early-type massive star and a compact companion, neutron star 
(NS) or black hole (BH), on a close orbit. Accretion of stellar wind 



Table 1. Stellar mode l parameters (see the full set of parameters in 
iFeldmeier eTal] il997al)) 



Spectral type 




09.7Ib 


Mass 


Af. 


34 Mq 


Radius 




24i?0 


Terminal speed 




1850km/s 


Mass loss rate 


M 


3 X 10"*^ A/0 yr 



^ IWaldron & Cassinellil )2010l) pointed out that the extreme UV radiation 
may be important for the mass-loss diagnostics based on resonant lines. 



onto the compact companion powers the high X-ray luminosity of 
~ 10^^ -10^^ erg s~\ HMXBs can be divided in subclasses de- 
pending on the spectral type of the donor star - supergiant (SG) or 
Be-star, and the type of variability - persistent or transient. 

In this paper we consider SG HMXBs. The basic physics of 
neutron-star a ccretion in stellar winds wa s described in the sem- 
inal paper by [Davidson & Ostrikeil ( Il973h . who suggest that the 
Bondi-H oyle-Lyttleton accretion po wers the X-ray luminosity in 
HMXBs. Ik^ man & McCravl ( fl982l) presented theoretical models 
for the temperature and ionization structure of spherically symmet- 
ric, constant density, gaseous nebulae surrounding compact X-ray 
sources. Two-dimensional hydrodynamic simulations of the gas 
flow in the o rbital plane of a mas sive X-ray binary system were 
presented bv iBlondin et al. I ( 119901) . The simulations revealed the 
presence of dense filaments of compressed gas formed in the non- 
steady accretion bow shock and wake of the compact object. The 
underlying assumption in all these models is an initially smooth 
donor stellar wind with the CAK density distribution. In a recent 
review, Edgar ( 2004}) demonstrated that despite the simplifications 
made, Bondi- Hoyle-Lyttleton pred ictions for the accretion rate are 
quite accurate. iNegueruelj j20IOl) reviewed recent works on stel- 
lar wind accretion, and suggested that the real sources differ from 
the idealized situation of Bondi-Hoyle-Lyttleton accretion in three 
main ways: the winds of massive stars are highly structured; accre- 
tion to a compact object is an unstable process; and the magnetic 
field of the neutron star may affect the flow of material. In this paper 
we investigate the immediate consequences of stellar wind clump- 
ing for the accretion rate, wind ionization, and the absorption of 
X-rays in the unperturbed portion of the wind. 

The observatio nal evidence of clu mped accretion flow in SG 
HMXBs is growing. ISako et al. I ( l2003l) reviewed spectroscopic re- 
sults obtained by X-ray observatories for several wind-fed HMXBs. 
They concluded that the observed spectra can be explained only as 
originating in a clumped stellar wind w here cool dense clumps are 
embedded in rarefied photoionized gas. Ivan der Meer et af] j2005h 
studied the X-ray light curve and spectra of the HMXB 4U 1700- 
37. They showed that the feeding of the NS by a strongly clumped 
stell ar wind is consis tent with the observed stochastic variability. 

lin't Zaiidl | |2005|) proposed that the flaring behavior of super- 
giant fast X-ray transients (SFXT) can be explained by the accre- 
tion of wind blobs. This suggestion was corroborated by the stud- 
ies ofhard_2C;ray_^ares_and^ emission of SFXT systems 
bv IWalter & Zurita HerasI Emh. who estimated wind clump pa- 
rameters and found them to be in agreement with the ma croclump- 
ing paradigm of stellar winds. iNegueruela et alj 1 I2OO8I) presented 
a common framework for wind accreting sources, in the context 
of clumpy wind models. They postulated that the different flaring 
behavior can be an immediate consequence of diverse orbital ge- 
ome tries. 

iKrevkenbohm et al] ( l2008h studied the temporal variability of 
Vela X- 1. They interpreted flares and off states as being due to a 
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Figure 1. A snapshot of the density (upp er panel) and the veloc ity (lower panel) structure of the stellar wind in a 09.5 supergiant star as predicted by 
time-dependent hydrodynamic simulations jpeldmeier et al.|[l997al lbl). The upper panel shows that the density fluctuations in the wind have many orders of 
magnitude, while the lower panel shows the strong wind velocity jumps. 

strongly structured wind of the optical companion. They suggested 
that when the NS in Vela X-1 encounters a cavity with strongly re- 
duced density, the X-ray flux drops, thus triggering the onset of the 
propellbibdefinitions.texer effect, which inhibits further accretion 
( off states). From t he analysis of the flaring behavior of VelaX- 
l lFurst etalj tOld) concluded that a mixture of a clumpy wind, 
shoc ks, and turbulence can explain the flare brightness distribu- 
tion. ISidoli et all ( 1201 Ol) suggested that the X-ray variability of the 
H MXB IGR 108408^5 03 may be explained by stellar wind clump- 
ing. |Bozzo|etal](|2oTl|) observed a bright X-ray flare in the HMXB 
IGR J18410-0535. They concluded that the flare was produced by 
the accretion of massive clump onto the compact object hosted in 
this system. 

iDucci et alj ( |2009|) developed a phenomenological stellar 
wind model for OB supergiants to investigate the effects of accre- 
tion from a clumpy wind on the luminosity and variability prop- 
erties of high-mass X-ray binaries. They concluded that the model 
can reproduce the observed light curves well. They also pointed 
out that the overall mass loss from the supergiant should be about 
a factor of 3 - 10 smaller than the values inferred from ultraviolet 
line studies that assume a homogeneous wind. This, as well as all 
other works on the phenomenology of clumped winds in HMXBs, 
account only for the presence of dense wind clumps and ignore the 
non-monotonic wind velocity. 

In this paper we employ time-dependent hydrodynamic sim- 
ulations of non-stationary stellar winds and use the predicted non- 
monotonic density and velocity fields to compute the accretion rate. 
We deliberately concentrate on the general consequences of stel- 
lar wind clumping for the HMXBs. Multi-D hydrodynamic models 
and the processes in the immediate vicinity of the compact object 
are beyond the scope of this paper We consider here systems of in- 
termediate and small X-ray luminosity, such that the donor stellar 
wind at large is not perturbed by the presence of the compact object 
and its X-ray emission. 

The paper is organized as follows. In Sect.|2]we describe the 
hydrodynamic models used to obtain the time-dependent wind den- 



sity and velocity. The Bondi-Hoyle-Lyttleton accretion formalism 
is applied in Sect.[3]to simulate the resulting time-dependent X-ray 
luminosity. The photoionization in the clumped wind is considered 
in Sect.|4l In Sect.|5]we apply our 2D model of a clumped stellar 
wind to predict the changes in absorption during the orbital motion 
of the NS. The discussion is presented in Sect.|6]and conclusions 
are drawn in Sect.|7] 



2 HYDRODYNAMICAL SIMULATIONS 

We employ the results of hydro dynamical simulations o f the line 
driven stellar wind presented in iFeldmeier et al. I( ll997j|bl) . These 
hydrodynamic simulations derive the dynamic and thermal struc- 
ture of stellar winds from the basic underlying physical principle, 
namely the acceleration of stellar wind by the line scattering of 
stellar UV photons (CAK). 

The simula tions were performed using the smooth source 
function method i lOwocki& Puis! 19961) . The line-driven instability 
(LDI) is handled by a careful integration of the line-optical depth 
with a resolution of three observer-frame frequency points per line 
Doppler width and by accounting for non-local couplings in the 
non-monotonic velocity law of the wind. The angle integration in 
the radiative flux is limited to one single ray which hits the star 
at ~0.7-R, . This accounts with sufficient accuracy for the finite- 
disk correction factor. The pure hydrodynamics part of the code 
is a standard van Leer solver. The seed perturbations for unsta- 
ble growth, in the shape of turbulent variation of the velocity at 
a level of roughly one third of the sound speed, are introduced at 
the base of the wind. These perturbations have a coherence time 
in the friction term of 1.4 hr, which is not too far off the acous- 
tic cutoff period of the star at which acoustic perturbations of the 
photosphere sho u ld accumulate. Further details can be found in 
lFeldmeiejfl995l) J Feldmeier et al. I( ll997allbl) . 

The hydrodynamic simulations were performed using stel- 
lar parameters of a typical late O-type supergiant star COrionis 
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(09.7Ib) (see Table[T). A snapshot of a wind structure at some given 
moment of time is shown in Fig.[T] The variety of dynamical struc- 
tures in the instabihty-induced, line-driven wind turbulence seems 
to be similarly intricate to that found in other turbulent flows in as- 
trophysical MHD settings. There are indications of the presence of 
a quasi-continuous hierarchy of density and velocity structures in 
the wind. Dense shells are formed close to the star, in a first stage of 
unstable growth. These dense shells have rather small radial extent, 
and contain a large fraction of the wind mass. In the framework 
of these hydrodynamic models, the colloquial term "wind clumps" 
refers to these dense sh ells. The space between the shells is filled 
by tenuous and hot gas l lFeldmeier eTalll 1997allbh . 

Further out in the wind, in a second stage of unstable growth, 
a somewhat different type of dense structures are formed: the small 
"clouds" that are "ablated" from the outer rim of the remaining 
mass reservoir at CAK densities, and are accelerated by the stel- 
lar radiation field through the emptied regions, eventually colliding 
with the next-outer shell. The X-ray emission from single early type 
star winds is explained by these cloud-shell collisions. 

While the average velocity field follows the so-called /3-law: 



v{r) 



1/t) 



(1) 



where Voo is the wind terminal speed, there are strong velocity 
jumps with negative gradients across the shells. In Figure|2] we 
show the wind velocity at some specific radius in the wind as func- 
tion of time. The figure illustrates, that relatively close to the star 
(we selected a radius 2.5R, from the stellar surface) the wind ve- 
locity changes by up to an order of magnitude on a time scale of 
hours. 



3 ACCRETION OF A NON-STATIONARY STELLAR 
WIND 

3.1 Basics of Bondi-Hoyle accretion 

The fo llowing is based on the treatment of [Davidson & Ostriked 
( Il973h . They considered a neutron star of mass mx, traveling at 
relative speed iirci through a gas den sity p. If Vrei is super sonic, 
the gas flows in a manner described bv lBondi & Hovlel ( Il944h . The 
mass accretion rate is given by 



where 



' acci' 



2Gm> 



(2) 



(3) 



and the factor ( is intended to correct for radiation pressure and the 
finite cooling time of gas. In mo derately luminous X-ray sources 
C<1 ( iDavidson & Ostrikei|[T973h . For our calculations we adopt 

C = i- 

In non-stationary stellar winds, the wind density p(a, t) and 
the wind velocity v„{a,t) depend on time (see Section|2]l. Com- 
bining equations ^ and ^ the accretion rate is 
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-p[a,t). 



The relative velocity 



(4) 



(5) 



where is the stellar wind velocity at the neutron star orbital sep 
aration a, and the orbital velocity of the neutron star is given by 
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Figure 2. The model (see text) wind velocity as function of time at the 
radius a = 2.5 iJ* . The amplitude of velocity variations may exceed order 
of magnitude in only few hours. 



The X-ray luminosity of an accreting neutron star is then 



(7) 



2 



(6) 



where c is the speed of light, and is a constant depending on the 
exact physics of accretion and typically taken to be ry ~ 0.1. 



3.2 Syntiietic accretion light curves 

The time dependent hydrodynamic simulations described in sec- 
tion|2]provide absolute values for the wind density and velocity as 
functions of radius and time for a star with parameters as shown 
in Table[T] We use these density and velocity to compute synthetic 
X-ray light curves according to Eqs. (|4l-(|7](. 

Adopting a mass of the compact object mx = 1.4 Mq and an 
efficiency parameter = 0.1, we generated simulated light curves 
for different values of the orbital separation a. The resulting pre- 
dicted variations in X-ray luminosity Lx are shown in Fig.jS] 

The synthetic light curves are highly variable, with X-ray lu- 
minosity changing by up to eight orders of magnitude! This strong 
variability is a consequence of the highly variable density and ve- 
locity in the stellar wind. Specifically we would like to empha- 
size that the shocks and corresponding jumps in the velocity with 
Sv ~ 1000 km s~^ are responsible for the predicted strong vari- 
ability, since the accretion rate scales with oc v^"^. Thus models 
which account only for the density fluctuation s in the winds (e.g. 
IWalter & Zurita Herasl2007l : lDucci et al.l2009l) may underestimate 
the resulting variability. 

It is interesting to note, that the character of the light curves 
cannot be characterized as steady quiescent emission and some 
strong flares on top it. Reflecting the nature of non-stationary in- 
homogeneous stellar winds the accretion light curves are truly 
stochastic. 

The X-ray observations of real celestial HMXBs are sensi- 
tive to the X-ray flux only about some threshold. In an arbitrary 
way we choose this threshold corresponding to the X-ray luminos- 
ity 10'^'* erg s~^. Figure|4] shows the fraction of time the luminosity 
of model source is above this threshold. It demonstrates that the 
fraction of time the source is X-ray bright decreases sharply with 
increasing orbital separation. In HMXBs with larger orbital separa- 
tion, the periods of highest accretion are significantly less frequent. 

We calculate the average X-ray luminosity resulting from di- 
rect accretion, L^, the standard deviation ctl, and the relative vari- 
ability aL/Lx- Our models show that there is no apparent depen- 
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Figure 3. Synthetic X-ray light curves for Bondi-Hoyle accretion of a non-stationary stellar wind onto a neutron star with mx = 1.4 A/q. The supergiant 
donor parameters are given in TabIe[T] From top to bottom panel, the orbital separation a is decreasing, as indicated in the lower left corner of each panel. The 
average X-ray luminosity is indicated in the upper left corner. The relative variability / L'^ = 4.8, 4.9, 4.9, 4.6 from top to bottom panel. 
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Figure 4. The fraction of time the synthetic accretion hght curve is above 
the level 10^^ erg s^^ in dependence on the orbital separation a. 



dence of the relative variability on orbital separation within the pa- 
rameter space we have chosen. This reflects the stochastic nature of 
stellar wind even relatively far from the stellar core. 



4 THE PHOTOIONIZATION PARAMETER IN THE 
CLUMPED WIND OF HMXBS 

The total X-ray emission in HMXB results from direct accretion 
which we considered in the previous section as well as from the 
photoionized fraction of stellar wind. Due to the highly variable 
accretion rate, the ionizing source is highly variable (see Fig.|3]l. 
Moreover, the ambient stellar wind is highly inhomogeneous (see 
Fig.[TJ. In this section we investigate how the ionization parameter 
is changed in a clumped wind compared to a smooth medium - 
the latter is commonly adopted in the computations of spectra from 
photoionized plasma. 

To study these effects we adopt a statistical approach, assum- 
ing a large number of clumps. The use of a statistical approach is 
convenient for at least two reasons. First, the non- stationary hy- 
drodynamic simulations we presented in the previous section are 
ID, and while they are well suited to compute the luminosity of 
the accreting point source, we need additional information to study 
spatially extended photoionized regions. Seco nd, there are sophis- 
ticate d photoionization codes, such as xstar dKallman & Bautist3 
I2OOII) that account realistically for atomic physics of photoionized 
gases. Using a statistical treatment of the clumped stellar wind, we 
attempt to find correction factors that would allow to adjust the in- 
put parameters for a photoionization code in order to reflect the 
wind inhomogeneity. 

The fundamentals of interaction of X-ray sourc e s wit h 
ambient stellar winds are developed in Tarter et al. I ( Il969t) . 
iTarter & Salpeted ( fl969t) . lHatchett & McCravl ( Il977l) considered a 
model of a constant X-ray source immersed in a stellar wind from 
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its companion star. Tiie stellar wind model was as predicted by the 
then recent CAK theory. These works consider a nebula gas of con- 
stant density A'' (or smooth density gradient) illuminated by a con- 
stant source of X-ray radiation. The ratio N^/N is determine d in 
each case from charge neutrality. Adopting iTarter et al.l j 19691) no- 
tations, neglecting absorption one may take for the mean primary 
photon intensity at frequency v ata distance r from the source 



J{u)du ■ 



L 

47j-y.2 



(8) 



where function /{v, kTyi) describes the spectral shape of the X- 
ray radiation. The equation of ionization equilibrium for a homo- 
geneous medium is 

N{Z,z + l)N,a{T,) = N{Z,z) j J[v)a{u)Au + N,C{T,), 

(9) 

where C is the coUisional ionization rate and g is the photo- 
ionization cross-section for the changing ionic charge 2 — >■ (z + 1) 
of ion with nuclear charge Z\ a is the total coefficient for recombi- 
nation for (2 + 1) — z. 

As follows from Eq. ^ and Eq. ^ for an optically thin gas 
in local ionization and thermal balance, irradiated by a constant 
point source of X-rays of given spectral shape, only a single pho- 
toionization parameter must be specified in order to determine the 
ionization state and temperature of any local region of the gas: 
^ = Lx/Na^, where N is the local atomic number density of 
the gas at the orbit of the X-ray source, and a is the distance from 
the center of the primary to the X-ray source. 

In the medium where the density is inhomogeneous, Eq. ^ 
shall be rewritten as 

{N{Z, z+l)N,)a{n) ^ {N{Z,z)) J J{u)a{u)du+{N,)C{n), 

(10) 

where the averaging over some local volume is represented by the 
angular brackets. Considering Eq. l llOt it is clear, that, in general, 
the ionization parameter ^ is not the same as in smooth density 
case. It is also clear, that in the medium with density fluctuations, 
the number of recombinations which is proportional to the den- 
sity squared is enhanced compared to the number of ionizations 
which is proportional to the density. Therefore, in the clumped stel- 
lar wind, the photoionization parameter must be reduced compared 
to the smooth wind. Moreover, the size of the photoionized region 
will be reduced in a clumped wind compared to the smooth wind 

of the same mass-loss rate. 

To describe the density fluctuations we follow I AllenI h973l) 
and define the wind inhomogeneity parameter 



(11) 



Then the root mean square (rms) local density is A^rms = {N)\/X. 
The rms density describes in a statistical sense the magnitude of 
density variations. 

Assuming that there are no dramatic changes in the ionization 
state over some local volume, we can write 



{N{Z,z + l)N,) ^ X{N{Z,z + l)){N,), 



(12) 



In this case, the photoionization parameter for a clumped stellar 
wind is 



X{N)c 



(13) 



i.e. it is reduced by factor X compared to the smooth wind case. 

The particle density averaged over some unite volume at dis- 
tance r, (A'^) is defined by the continuity equation: 



(N) 



M 



(14) 



where /i is the mean molecular weight, ttih is the atomic mass unit, 
and v{r) is a stationary velocity law. 

The simplest structure of a clumped wind is a two phase 
medium where all matter is clumped with dense clumps occupying 
only a small volume compared to the interclump medium filling 
the large volume. The interclump medium can be assumed void or 
filled with tenuous gas. In this case a volume filling factor /v can be 
defined as the fraction of total volume Kot that is filled by clumps: 
/v = Ki/Vtot. If the interclump medium is void X = f-^^. The 
clump density in such two phase medium is A^c = X{N). If the 
interclump medium is not void, its particle number density can be 
approximated as: 



Ni,^{l-XfY){N). 



(15) 



4.1 On the clump optical depth 

Equation ([8) neglects the absorption of ionizing radiation. This con- 
dition may not be valid when clumps are optically thick for the 
ionizing radiation. In this section we investigate the restrictions on 
the clump parameters imposed by this condition. To do so we es- 
timate the geometrical size a clump which has optical depth unity. 
Assume for simplicity that the interclump medium is void, so that 
/v ~ X^^, the optical depth of an isotropic clump (i.e., with the 
same dimensions dc in 3D) for the incident radiation at wavelength 
A is 



Tc{\) = p^Kxfv^dcR, 



(16) 



where R, [cm] is the stellar radius, dc is in the units of i?, and 
K,\ [cm^ g^^] is the mass absorption coefficient. The density of the 
cool wind (pw) can be defined from the continuity equation. 

The clump dimension dc can be derived from Eq. il6l . Adopt- 
ing rc(A) = 1: 



/v 



r.(A) 



(17) 



where r, is a parameter and a /3-law for the velocity is assumed. If 
at some radius r, the clump size is larger than d^i^^p, such a clump 
is no longer optically thin for X-ray radiation at A. The parameter 
r, (A) can be conveniently expressed as 



r.(A) ^ 722 



M- 



(18) 



where M-q is the mass-loss rate in units 10 ® Mq yr ^, is the 
termi nal velocity in [kms~^], and 7?., — Rt/RQ jOskinova et al.l 
l20Ilh . 

The K-shell opacity of stellar wind is a strong function of 
wavelength, and can be taken as ~ «;o-^^> with 7 in the range 
2-3, and kq a constant that de pends on abundances dSaum et aP 
Il992l : llgnace & OskinovJl999l) . It follows from Eqs. GT} and (TDl 
that the same clump can be optically thin at shorter wavelengths 
and optically thick at longer wavelengths. 

Realistic values of the mass-absorption coefficient have to be 
computed usi ng non-LTE model stellar a tmosphere codes, such as. 



e.g., PoWR Ih amann & Grafenen 



sample of O-stars are given in e.g 



|2004). The va l ues of kx for a 



Oskinova etallfeOOeh 
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The SG HMXB are hard X-ray sources. The maximum 
of spectral energy dis tribution in VelaX-1 is at ~ 4 A (e.g. 
IPoroshenko et al.ll201 lb . Using parameters of Ori and assuming 
/v = 0.1, at a = 2 7?., dJi=^p(COri) ^ 20«:^\ Le. at the dis- 
tance 2 R, from the stellar core, only geometrically large clumps 
with diameter larger than ~ 0.5 7?* are optically thick for the ion- 
izing radiation shortwards of 10 A. 

Such a clump size appears rather large, therefore it is probably 
reasonable to assume that the wind clumps are optically thin for the 
ionizing radiation resulting from the accretion onto a NS0 



4.2 Optically thin clumps 

As shown above, the wind clumps are most probably thin for the 
hard X-ray emission resulting from the accretion on a compact ob- 
ject. In this case, the "microclumping" approximation can be used. 
If the mean primary photon intensity is described by Eq. ((Sj, then 
in a two phase medium where X — /y^, the photoionization pa- 
rameter can be written as 



{N)c 



(19) 



The meaning of Eq. l |I9| l is that the ionization parameter in a 
clumped wind is reduced by a factor /v whe n the average density 
at dist ance a is given by Eq. l ll4t . Following iHatchett & McCravl 
( Il977b the nondimensional variable q = £,{N)a'^{Lxfv) ^ de- 
scribes the shape of the surfaces for which ^ = constant. In 
clumped stellar wind, the parameter q will be larger than in a 
smooth wind. Considering the case of constant wind velocity, it 
would mean that the radius of spherical photoionized region is 
smaller in clumped wind compared to the smooth wind. This is 
easy to understand recalling that the number of recombination is 
enhanced in clumped wind compared to smooth wind, while the 
number of ionizations remains the same. 

Over the last decade there have been numerous stu dies to es- 
tablish the degree of wind clumping in OB supergiants. IPuIs et al] 
( l2006h analyzed a sample of 19 Galactic 0-type supergiants/giants, 
covering spectral types 03 to 09.5 using multiwavelength observa- 
tions in radio, IR, and optical. They used the conventional assump- 
tion of the interclump medium being void. It was found that the 
upper limit on the volume filling factor in the innermost region ( 
r < 2 7?*) of stars with in emission is about 0.25. It is likely 
that the winds of donor stars in HMXBs have similar properties. 
In this case, the ionization parameter should be reduced by a mini- 
mu m factor of four to acco unt for wind clumping. 

iFullerton et Zl l l2006h analyzed 40 Galactic O-type stars by 
fitting stellar wind profiles to observations of the Pv resonance 
doublet in the UV. The found that to reconcile mass-loss rates em- 
pirically infened from the analysis of UV and radio/optical spec- 
tra, the volume filling factors must be /v <0.01. Similar conclu- 
sions were reached bv lBouret et al.l ( l2005b who found /v ~ 0.04 
in an O-type supergiant and fy ~ 0.02 in an 0-type dwarf. If 
these clumping factors were realistic, the ionizati on parameter ^ 
would have to be reduced by factors of 25 ... 100. lOskinova et al.l 

^ Note, that the intrinsic X-ray radiation in single SG OB stars is thermal 
and soft, with the maximum of t he differential emissi on measure distribu- 
tion at s:^ 60 A or even softer (e.g. lRaassen et al.l2008l) . Even geometrically 
very small clumps with diameter io~3-..— 5 ^re optically thick at these 
wavelengths. Therefore, generally, in the inner wind regions of single SG 
OB stars the optically thin clump approximation for the intrinsic X-ray ra- 
diation is false. 



( l2007b showed that correctly accounting for the wind porosity in 
the analysis of UV resonant lines eliminates the need for very 
small volu me filling factors and t hat realistic values are rather 
/v ~ 0.1. Sundqvist et al. 1 201 1 ) also used the hydrodynamic 



models of IFeldmeier et al. ( 1997al lbl) and found that in these mod- 



els the typical clumping factors are 



10. In this case, in the 



studies of the X-ray spectra from photoionized gas, the ionization 
parameter ^ should be reduced by a factor <10 in clumped wind 
models compared to the unclumped ones. 



4.3 Clumps of arbitrary optical depth 

While it is likely that the wind clumps are optically thin for the 
hard ionizing radiation, (see Section l4~]l . for completeness in this 
section we consider a case of optically thick clumps. Accounting 
for the clumps of arbitrary optical depth in stellar atmosphere mod- 
els is sometimes referred to as "macroclumping". If the optically 
thick structures, Tc^JI, are present in the wind, the mean primary 
photon intensity of ionizing radiation must be corrected for the 
absorption in the clumps 



J. 



fikT) . 



f(r) 



(20) 



where f{kT) is a function of atomic parameters and tempera- 
ture that describes the spectrum of the incident ionizing radia- 
tion, and r is the effe ctive optical depth of the clumped medium 
( IFeldmeier et alj2003l) . 

At a distance 7?x from the source of ionizing radiation the 
effective optical depth is 



nccrc(l - e '^°)drx, 



(21) 



where nc{r) oc v^^{r)r^'^ is the average number of clumps 
within a unit volume, and ac is the clump geometrical cross- 
section jOskinova et al.ll2008ll201 lb . Note that the stellar wind pa- 
rameters M and k\ appear in the expression for r only indirectly 
via Tc given by Eq. il6\ . Besides stellar wind parameters, the ef- 
fective optical depth also depends on the geometrical distribution 
of the clumps via nc and ctq- The clump cross-section is different 
for the different clump geometries, e.g. for the spherical clumps 
(Jc oc dQ, while for the oblate clumps the cross-section depends on 
the clump orientation. In the limit of rc S> 1, the effective optical 
depth is fully determined by the geometrical distribution of clumps 
( IFeldmeier et alj2003l : lc)skinova et all2004b . 

Now we would like to address the question of ionization bal- 
ance in the wind where the largest portion of the matter is in form 
of optically thick clumps. The X-ray emission from such medium 
will have two components: X-rays due to scattering of radiation 
within the clumps, and the X-rays from the photoionized inter- 
clump mediu m. The scattering of pho tons in such a medium was 
considered bv lSuleimanov et aU ( l2003h . Here, we consider the ion- 
ization of the interclump medium and neglect the diffuse radiation. 
Using Eq. ijTSj for the particle number density in the interclump 
medium, the ionization parameter for the tenuous interclump gas 
becomes 



(22) 



{iV)a2 1 - Xfv ■ 

When the bulk of stellar wind mass is in the form of a few dense 
optically thick clumps, the interclump medium can be very tenuous 
X — fy^.ln this case the effective optical depth can be very small 
(i.e. the wind is very porous). As follows from Eq. l l22t . in this case 
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Figure 5. Sketch of the coordinate system for the fragmented wind model. 
X-ray radiation with intensity is emitted at the location with coordinates 
(pi, The cool wind consists of clumps that are randomly distributed. 
For clarity, only one clump (shaded area) is shown in the sketch. The X-ray 
radiation propagating towards the observer at 2 — > oo has to cross a clump. 
The optical depth across the clump for the ray in z direction is Tc , reducing 
the line intensity /J^ by a factor exp(— Tc). 



the ionization parameter in the tenuous interclump gas can be very 
high, much higher than would be expected for the smooth wind of 
the same mass-loss rate. 

We conclude that clumping strongly affects the ionization bal- 
ance. A higher degree of ionization in the less dense interclump 
gas within the photoionized region can be expected. The size of the 
photoionized region will be much larger in case of optically thick 
clumps compared to a smooth wind. 




" 5 » V V 



Figure 6. One random trial of our 2-D model where the clumps are ran- 
domly distributed in radius (uniformly between 1.2 and 316 On av- 
erage there are 400 clumps along each radial direction (see text for de- 
tails). Only the inner pai't of the wind between li?, and 5i?« is shown. The 
smooth wind is in-between 1 i?,, and 1 .2/? * . The filled circle represents the 
stellar core; the outer circles are at 2i?* and AR, to scale. This random 
model was used to calculate optical depths shown in Fig.[8]and Fig. I?] 



5 ABSORPTION OF X-RAY RADIATION IN A 
CLUMPED STELLAR WIND 

The X-ray photons emitted either in close vicinity to the NS, or far- 
ther out in the photoionized portion of the stellar wind, suffer the 
K-shell continuum absorption in the cool component of the stel- 
lar wind as they propagate towards the observer. As discussed in 
Sect. l4.1l the wind clumps are, most likely, optically thin for the 
radiation at the short wavelengths (shortward of ~ 8 A). On the 
other hand, the wind clumps or larger scale wind structures can be 
optically thick for the softer X-rays. These optically thick inhomo- 
geneities can be traced observationally as changes in the absorb- 
ing column density on the time scale comparable to the orbital pe- 
riod. Such variab ility of absorbing column is indeed detected (e.g. 
iFiirst et alj201 ih . If the stellar wind of donor star were smooth, the 
absorbing colunm would smoothly increase from the moment when 
the NS is at the periastron to when it is in the apastron. 

In this section we investigate the K-shell continuous X-ray 
absorption in stellar wind. Our goal is to model the absorption 
of monochromatic X-ray flux emerging from the NS in depen- 
dence on orbital p hase. We employ a 2 D model of a stochas- 
tic stellar wind by lOskinova et alj ( |2004|) . where the model de- 
tails are fully described. Here we only briefly introduce its ba- 
sic features. The model assumes the structure and physical con- 
ditions in the w i nd acco rding to the hydrodynamic simulations by 
iFeldmeier et al. Hl997al lbl). The emission is decoupled from absorp- 
tion, which allows for the formal solution of the radiative trans- 
fer equation. In order to take advantage of the rotational sym- 
metry of our model around the z axis, cylindrical coordinates 
(p, z, tfi) are used. Along each line of sight towards an X-ray emit- 
ting region located at spherical coordinates (ri,i?i), the emergent 
monochromatic X-ray intensity I^{p, z, Lp) is reduced by a factor 



exp(-r(p,2:,v?)): I+{p,z,ip) = /"(p, z, ^) exp (-r(p, lys)). 
The optical depth r(p, z, ^p) is given by the summation of optical 
depths of all clumps on the line of sign. The geometry of the model 
is shown in Fig.[5] Due to the lack of synchronization between dif- 
ferent radial directions, the fragments probably have rather small 
lateral extent. Despite the expansion due to the internal pressure, 
the cool fragments are maintained to distances of hundreds of the 
stellar radius, before they gradually dissolve into a homogeneous 
outflow. Our 2-D stochastic wind model is designed to describe in 
a generic way the structures and physical conditions of the inho- 
mogeneous stellar wind: 

1. the flow is spherically symmetric in a statistical sense and prop- 
agates according to a /3 velocity law. The mass-flux and total radial 
optical depth in each radial direction is preserved. 

2. all clumps cover the same solid angle as seen from the center of 
the star, and preserve this angle during their propagation. 

3. we consider two different shapes of clumps, the spherically sym- 
metric clumps ("balls") and radially compressed clumps ("pan- 
cakes"). 

4. to simulate stochastic X-ray emission we model X-rays as dis- 
crete random flashes of X-ray light along the NS orbit. There is no 
self-absorption by emitting material and no re-emission of X-rays 
which became absorbed by cool fragments. 

5. we consider an edge on system, i.e. the observer is located in the 
orbital plane of the NS 

We neglect the density fluctuations which are likely to result in 
the NS wake. The hydrodynamic simulations of wind structure that 
predict the existence of a wake were performed for smooth winds 
dBlondin et alj 19901) . It is not not obvious how accounting for wind 
clumping can affect the results of hydrodynamic simulations and 
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Figure 7. The variation of the factor exp {—t\) in dependence on orbital 
phase for the X-ray radiation at 8 A (1 .55 IceV) emitted by a source located 
at 2 . The horizontal axis shows the orbital location of the X-ray 
source, i?;. At 0° the source is in conjunction, while at 180° it is totally 
eclipsed. The red curve shows the changes of cxp (— t;^) in a smooth wind. 
The parameters of the wind are given in Table[T] A velocity law with expo- 
nent /3 = 0.8 is assumed. The stochastic wind model is calculated using 
the random trial shown in Fig.[6] Upper panel: the wind model where the 
clumps are assumed to be spherically symmetric {'"balls"). Lower panel: 
the wind model where the clumps are assumed to be flattened in the radial 
direction {'"pancakes"). 



the shape and density in the wake. Moreover, the presence of the 
wake should introduce some regular pattern in the absorption of 
X-rays along the orbital motion of the NS, while our goal is to 
investigate the basic effects of stochastic winds on the absorption 
of X-rays. 

We also neglect the changes in the ionization structure of the 
wind due to the X-ray photoionization close to the NS. As dis- 
cussed in Section|4l the clumping is likely to reduce the size of 
photoionization zone. The optical depth of the clumps within the 
photoionized portion of the wind may be different from the rest 
of the wind, however this effect may be small on average due to 
the difference in the size of the general wind, and its photoionized 
portion. 

Figure|6]shows a snapshot of the inner wind structure, between 
1 i?, and 5 i?, (in the simulations clumping is maintained up to 
316 -R*). We assume the lateral extend of clumps to be 1°. To study 
how the wind optical depth changes during the orbital motion of the 
NS, we assumed the donor stellar parameters as listed in Table[T] 
The clump number density obeys the continuity equation, therefore 
the clump density is higher in the innermost part of the wind and 
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Figure 8. The same as in Fig .[T] but for the radiation emitted at 19 A 
(0.65 keV). 



gradually decreases outwards. To find the optical depths towards 
the observer at different orbital location, we sum up the contribu- 
tions from all the clumps along the line of sign. 

The mass-flux in each direction is preserved in our model. To 
calculate the optical depth of each clump (in total there are 144 000 
clumps in the wind model, part of which is shown in Fig.|6]l we 
assume that a clump mass equals to the mass of the smooth wind 
con fined between to random subsequent radii in the wind ra and rb 
(see lOskinova et all ( l2004h for details). For a star with the param- 
eters shown in Table 1 and solar abundances, the mass-absorption 
coefficient ka is ~ 40 cm^ g~^ at A = 8 A and ~ 170 cm^ g~^ at 
A = 19 A l lOskinova et"ai]|2006l) . 

Little is known about the geometrical shape of the clumps. The 
comparison with the observed X-ray emission line profiles in spec- 
tra of single O-stars indicate that th e clumps can be flattened in ra- 
dial direction jOskinova et al .120061) . The shape of clumps plays an 
important role in the transfer of X-ray radiation through the wind: 
in case of flattened clumps the wind opacity becomes anisotropic, 
while it remains isotropic (the same as for the smooth wind opac ity 
albeit smaller) in case of ball-like clumps dOskinova et"ai]|2006l) . 

Figures|7]and[8]show the transmission factor exp {—t\) for ra- 
diation at 8 A and 1 9 A ( 1 .55 keV and 0.65 keV, respectively) calcu- 
lated assuming different clump geometries and a smooth wind. The 
wind is more transparent for the radiation at shorter wavelengths. 
In case of flattened clumps the changes in optical depth are larger, 
because the clump optical depth depends on the clumps orientation. 

The changes in the wind absorbing column at different orbital 
phases should result in changes in the X-ray spectra. Indeed, the 
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Figure 9. The same as in Fig.[8] Upper panel: comparing the cxp (~t\) 
factor at 19 A (red) and 8 A (blue) at different orbital phases for a model 
wind that consists of flattened clumps. Lower panel: the same as upper 
panel but for a smooth wind model. 



strong variations of absorbing column density at diffe rent orbital 
phases are often deduced from the observed spectra (e.g lFiirst et al.l 
l20ld : lBozzoetalj201lh . 

The K-shell absorption in stellar wind is strongly wavelength 
dependent, changing by a two orders of magnitude across the 0.1- 
12 keV X-ray band. This is reflected in smaller variations in the op- 
tical depths predicted for the 8 A radiation (where the wind opacity 
is smaller, Fig.|7]l compared to the large optical depth variations at 
19 A (where the wind opacity is larger, Fig.[8]l. The attenuation of 
X-ray radiation in a clumped wind h as a weaker dependen ce on 
the wavelength than in a smooth wind jOskinova et al.ll2004l) . This 
prediction holds also for the HMXB. In Fig.[9]we compare the du- 
ration of the X-ray eclipse at two different wavelengths for clumped 
and smooth wind models. In a smooth wind, the eclipse at longer 
wavelength starts earlier, and lasts longer, while at shorter wave- 
length the eclipse is more narrow. On the other hand, in a clumped 
wind model, the difference between eclipse durations is less pro- 
nounced. Even in the phases close to the eclipse, optical depths 
may occasionally be as low as in the phases far from the eclipse. 



6 DISCUSSION 

It emerged from our analysis that the stellar wind clumping has 
drastic consequences for the phenomenology of SG HMXBs. The 
effect of wind structure on the accretion rate is profound. The hy- 



drodynamic simulations predict wind structure consisting of strong 
velocity and density jumps. Since the Bondi-Hoyle-Lyttleton ac- 
cretion rate depends on both velocity and density fluctuations, it 
strongly varies with time. Combining the Bondi-Hoyle model of 
accretion with hydrodynamic wind simulations we computed syn- 
thetic X-ray light curves for the different orbital separations (see 
Fig.©. 

Averaged over time model X-ray luminosity is in general 
agreement with observations. The moderate reduction of M by fac- 
tor 2 ... 3 compared to the unclumped wind models explains well the 
optical and the UV spectra of single stars, as well as average X-ray 
luminosity resulting from wind accretion in HMXBs. 

The application of the Bondi-Hoyle-Lyttleton approximation 
for the accretion rate predicts very strong X-ray variability: the X- 
ray luminosity changes by many orders of magnitude on the time 
scale of hours as a direct consequence of the strong density and ve- 
locity jumps in the accretion flow. The effect of the non-stationary 
wind velocity on the accretion rate was overlooked previously. Ne- 
glecting this effect may lead to erroneous estimates of clump mass 
and size from the total flux and duration of the X-ray flare. We em- 
phasize that the strong velocity jumps in the stellar winds lead to 
strong variability of accretion X-ray luminosity. 

The detailed comparison of model light-curves and the obser- 
vations will be a focus of our forthcoming study. However, it ap- 
pears safe to conclude, that the synthetic light-curves do not agree 
well with the observations over-predicting the amplitude of vari- 
ability. Given our simple approximation for the accretion this is not 
surprising. The detailed hydrodynamic studies that wind accretion 
is a highly complex process. E.g. when there is a density and ve- 
locity gradient in the flow, a transfer of angular momentum can be- 
come possible. The hydrodynamics of Bondi-Hoyle-Lyttleton ac- 
cretion in a mediu m with continuous de nsity and velocity gradient s 
was considered bv lRuffertI l fT997lll999h and lFoglizzo et all ( l2005h . 
To our knowledge, up to now there were no studies of accretions 
from a medium with strong velocity and density jumps. We spec- 
ulate that the detailed physics of the accretion flows regulates the 
temporal behavior of accretion and is pivotal in damping the strong 
variations due to the structured nature of stellar winds. 

Another process that can reduce the variability is the destruc- 
tion of clumps in the vic inity of the N S. As an interesting exam- 
ple of clump destruction. IPittardI ( |2007|) performed a hydrodynam- 
ical simulation of wind-wind collision in massive non-degenerate 
binary star. The simulations revealed that the clumps are rapidly 
destroyed after passing through the confining shocks of the wind- 
wind collision region (assuming adiabatic flow in the wind colli- 
sion region). Despite large density and temperature fluctuations in 
the postshock gas, the overall effect of the interaction is to smooth 
the existing structure in the winds. One may speculate that in con- 
ceptually similar fashion, the wind clumps and the wind velocity 
gradients are "smoothed out" in the strong shock of the accretion 
wak e in HMXB s^ 

iFiirst et al] ( 1201 Oh found that the accreted clump masses de- 
rived from the INTEGRAL data o n Vela X- 1 are on the order of 5 x 
10" - 10^1 g. |Bozzoetal.l l l201lh estimate a mass and a radius of a 
clump responsible for a flare in a SFXT as ~ 10^^ g and ~ 10^^ cm 
correspondingly. These are important independent estimates of a 
clump mass that can be compared with what we know about clumps 
in single stars. [Lepine & Moffat ( 1999) found stochastic line profile 
variations in the form of narrow (~ 1 00 km s ~ ^ ) emission sub-peaks 
on top of emission lines in Wolf-Rayet star winds. They explained 
their data by V^^NcliX^'^ "blobs" being presen t at any time in 
the line emission region. In lOskinova et alj ( I2OO6I) we reproduced 
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the Chandra observation s of 0-star X-ray Une pr ofiles by adopting 
a dynamic picture from iFeldmeier et alj ( Il997al lbl) hydrodynamic 
simulations - in each radial direction a clump is ejected once per 
dynamic time scale ta = R,/voo, which corresponds to a radial 
separation of one stellar radius when Voo is reached. Assuming that 
the wind consists of A^'c clumps, an average clump mass Mc is 

Mc ~ 4.4 X 10'^ ^^^^ [g] (23) 

Adopting a mass-loss rate M_6 ~ 1 [lO^^Af© yr~^], Voo ~ 
2000 [km s"^], and R, = 10i?Q with iVc ~ 10^ ... 10^ clumps 
in the wind, the average mass of a clump is Mc ~ 10^* ... 10^" g. 
Th ese number s comp are very well with those independently found 
by iFurst et al.l ( 1201 Ol) from the analysis of the INTEGRAL light- 
curves. 

In this study we briefly addressed the absorption of X-rays in 
clumped stellar wind. We found that the wind optical depth strongly 
varies in dependence on orbital phase. We neglect the photoioniza- 
tion wake. Partly, this can be justified, because not much is known 
about the driving of X-ray photoionized clumped winds. Classi- 
cally, the works that predict the decrease of driving force in the 
X-ray photoionized part of the outflow consider a smooth win d 
jMacGregor & Vitelloll983 : IStevenj|l988l : IFeldmeier et"ai]|l996h . 
In clumped wind the ionization balance may be altered, e.g. the 
low ionization stage that can drive stellar wind may be still present 
within the clumps. Our investigation of of changing optical depth in 
the wind around the compact object orbit demonstrates the poten- 
tial of high-quality X-ray observations of SO HMXBs to provide 
much needed detailed information about the parameters of wind 
clumping. E.g. the orientation of the clumps, their distribution, can 
be studied from the changing absorbing column in stellar wind, es- 
pecially during the X-ray eclipse. 

The important result of our study, is that the ionization param- 
eter is affected by the wind clumping. This follows from the basic 
fact that in the medium with density fluctuations, the recombina- 
tion process (p^ -dependent) have priority over the photoionization 
process (p-dependent). The clumping of stellar wind is an estab- 
lished fact, and the ionization parameter must be corrected for the 
clumping in the analyses of spectra of photoionized gas in HMXBs. 
The ionization structure, the temperature, and the size of photoion- 
ized region is different in the clumped stellar wind compared to a 
smooth wind with the same mass-loss rate. 



7 CONCLUSIONS 

1 . The results of time dependent hydrodynamical models of stellar 
winds predict a highly variable rate of wind accretion onto a com- 
pact object. The variability is largely due to the strong variations in 
the wind velocity (shocks). The highly non-uniform wind density 
also contributes to the variability of accretion rate. 

2. Combining the Bondi-Hoyle model of accretion with hydrody- 
namic wind simulations predicts variability on stronger levels than 
commonly observed. 

3. The X-ray light curves synthesized under the assumption of di- 
rect accretion predict strong variability of few orders of magnitude. 
The relative variability (jl/L'^ has no strong dependence on the 
binary orbital separation. 

4. The photoionization parameter depends on the properties of wind 
clumping. In optically thin case, the photoionization parameter has 



to be reduced by the factor X compared to the smooth wind model, 
where X is the wind inhomogeneity parameter. 

5. The absorption of X-rays in a clumped stellar wind strongly 
varies with the orbital phase. The degree of variability depends on 
the shape of the wind clumps, and is stronger for radially com- 
pressed clumps. 

6. In a clumped wind, the wavelength dependence of the duration 
of the soft X-ray eclipse is smaller compared to the smooth wind 
models. 

Overall, we conclude that the clumped non-stationary nature 
of the wind from the donor star strongly affects the X-ray light 
curve and leads to strong variability. 
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